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1. Introduction 

The nucleon-deuteron charge-exchange reaction is the subject of the investigation in 
the set of the experiments, which are started in the VBLHE JINR at STRELA and 
DELTA SIGMA |T] setups and in COSY j2] at ANKE spectrometer. These experiments 
are performed in special kinematics, where transfer momentum from the initial nucleon 
to the outgoing fast nucleon is close to zero. The goal of these experiments is to extract 
the additional information about the spin-dependent part of the elementary np — > pn 
amplitude from the nucleon-deuteron reaction. First this idea was suggested by Pomer- 
anchuk [3] already in 1951. Later, it was shown that in the plane-wave impulse approx- 
imation (PWIA) the differential cross section and tensor analyzing power T 2 q in the dp 
charge-exchange reaction are actually fully determined by the spin-dependent part of the 
elementary np — > pn amplitudes [2], [I]- 

However, the relative momentum of the two slow nucleons is very small under such 
kinematical conditions, where momentum of the emitted fast nucleon has the same di- 
rection and magnitude as the beam (in the deuteron rest frame). As a consequence, the 
final-state interaction (FSI) effects have to play very important role. How much is the 
FSI influence? It is the general question of our consideration in this paper. 

We consider the nd — > pnn reaction in kinematics of the DELTA SIGMA experiment, 
where the outgoing proton has the same direction as the projectile neutron and transfer 
momentum is close to zero. The neutron-beam kinetic energy changes from 0.8 up to 
1.3 GeV. The theoretical approach is based on the Alt-Grassberger-Sandhas formulation 
of the multiple-scattering theory for the three-nucleon system. We apply the matrix 
inversion method to describe the FSI contributions. Earlier, this formalism has been 
employed for description of the deuteron-proton breakup in the GeV region [5]- [6]. All 
calculations have been performed in the deuteron rest frame. The results, which obtained 
in the present paper, are very interesting, from our point of view, and useful for further 
investigations of this problem. 

The paper is organized as follows. The theoretical formalism is given in the Section 2 in 
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which we consider both the general description of the nd break-up reaction and the special 
case, where transfer momentum is close to zero. The results are presented in Section 3. 
The ratio of the nd — > pnn differential cross section to the free np scattering differential 
cross section at 6^ = has been studied. It is shown that final-state interaction is very 
important, although the FSI contribution decreases with the energy increasing. In Section 
4 the conclusions are given. 

2. Theoretical formalism 

In accordance with the Alt-Grassberger-Sandhas (AGS) formalism of the three-body 
collision theory [7]- [H], the amplitude of the neutron-deuteron charge-exchange reaction, 



U nd ^ pnn = V2< 123| [1 - (1, 2) - (1, 3)]£/oi|l(23) >= 5{p - Pl - p 2 - p 3 )J. (2) 



The state 1 1 (23) > corresponds to the configuration of nucleons 2 and 3 forming the 
deuteron state and nucleon 1 being the projectile, whereas the state < 123 1 represents 
the free motion of three nucleons after the reaction. The permutation operators for two 
nucleons (i, j) appear in this expression due to particle identity both in initial and final 
states. 

Iterating the AGS equations up to the second-order terms, we can obtain the expression 
for the transition operator Uqi. 



By definition, t 3 = t\ 2 is a two-nucleon transition operator (the others can be obtained 
via cyclic permutations). Here, go is the free nucleon propagator, 



n(p) + (2(0) -> p(pi) + n(p 2 ) + n(p 3 ) , 




is defined by the matrix element of the transition operator Uqi: 



Uqi = g 1 + h + h + t 1 g Q t 2 + hgoh + t 2 g t 3 + t 3 g t 2 + 0(t 3 ). 



(3) 



g = [E + zO - K] 
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with the kinetic-energy operator for the three-body system K and the on-shell energy 
E = Ei + E 2 + E 3 . Obviously, the term g^ 1 does not contribute to the on-shell amplitude 

Since we consider only such kinematical conditions, when outgoing proton carries away 
significant part of the initial neutron momentum, we can neglect the corrections due to 
the recoil reaction mechanism, where final fast proton leaves the deuteron without direct 
knock-out. In addition, we ignore terms, which correspond to the double scattering with 
participation of a fast nucleon. Then, the matrix element U n d^ pn n can be approximated 
by 

U nd ^ pnn = y/2<123\[l-(2,3)][l+t 23 (E-E 1 )g 23 (E-E 1 )]t s & m \l(23)>, (4) 

where the operator g 23 {E — Ex) is a free propagator for the (23)-subsystem and the 
scattering operator t 23 (E — Ex) satisfies the Lippmann-Schwinger (LS) equation with 
two-body force operator V 23 as driving term 

t 23 {E - Ex) = V 23 + V 23 g 23 {E - Ex)t 23 {E - Ex). (5) 

The operator tY 2 m is symmetrized NN-operator, tY 2 m = [1 — (l,2)]t 12 . 

The expression (J3j) can be schematically presented by the set of graphs (Fig.l). Here, 
two first terms correspond to the PWIA and the others represent the scattering with 
FSI of the slow nucleons. Let us rewrite the matrix element (jl]) indicating explicitly the 
particle quantum numbers, 

U n d^pnn = < Pi^Ti , p 2 m 2 T 2 , p 3 m 3 T 3 \ [1 - (2, 3)]u 23 tf 2 m \ pmr, l/lXM d 00 (23) >, (6) 

where u 23 = [1 + t 23 (E — Ex)g 23 (E — Ex)] and the spin and isospin projections are denoted 
as m and r, respectively. 

In momentum representation the deuteron wave function (DWF) tfjxM d (k) with spin 
projection is written as 

L 

|V> 1Md (k) >= £ J2 <LM L lM s \lM d >u L (k)Y^Ck)\lM s >, (7) 

L=0,2 M L =-L 
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with the spherical harmonics Yjf L (k) and the Clebsh-Gordon coefficients in the standard 
form. In our calculations, we use the pole parameterization of the deuteron wave function 



af+p 2 V 7T * A 2 + P 2 

Inserting in Eq.([2]) the unity 



1 = J dp'lp'm'r' >< p'm'r" 



and using the definition of the deuteron wave function, Eq.flZj), we get the following 
expression for the reaction amplitude: 

J = g < 2 m *2 ms \ SMs >< 2 m ' 2 2 m '^ SM ' s >< LM L lM s\ lM d >< \m"^m' 3 \lM s > x 



x Jdp 'lp ,SM s 



tS T '=ME-E x 
1 + m N - 



r> po', SM' S ML (|Po " q/2|)^ 3L (p^, - q/2) x 
Po — Po ~r *u / 

x < pimx, (po' + q/2) m 2 |t?To(£ - K) ~ *t=i(£ - ^)|pm, (p ' - q/2) m" > - 

- (2 - 3), (9) 

where is the nucleon mass, S corresponds to the spin of the two slow nucleons which 
participate in the final-state interaction. We have introduced here the momentum transfer 
q = p — pi, relative momenta p = |(pa — P3) and p ' = \{P2 — P3')- Henceforth, all 
summations over dummy discrete indices are implied. 

The two slow neutrons interaction is described by the wave function 

< ^poSm s tm t \Po'SM' s TM t >= 5{p - po')8 MsM > + 

+ 2 miV /2 : . n < PoSM s \t ST \p 'SM' s > . (10) 
Po 2 - Po 2 + «0 

This wave function contains the FSI part, which can be taken into account in different 
ways. 

In the present paper we use the matrix inversion method (MIM) suggested in [TQl [11] 
and applied to study the deuteron electro-disintegration [T2J [T5] and deuteron proton 
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breakup process [5j [6] . As in ref . [12] , we consider the truncated partial- wave expansion, 

< ^ P ~SM S TM T \P0 SM S™T >= $M S M' S S(P0 ~ Po') + (11) 

^max J 

+ E E Yl l (p )<lfiSM s \JMj>^(Po)<l'^SM' s \JMj>Y i r'(p ') 1 

J=0 Mj=-J 

where J max is the maximum value of the total angular momentum in nn partial waves 
and a = { J, S, T} is the set of conserved quantum numbers. 

Under kinematical conditions, where transfer momentum q = p — pi is close to zero, 
one can anticipate that the FSI in the 1 So-state is prevalent at comparatively small po- 
values. Then two-neutron wave function is 

< >= *(Po - Po') + ^Oo')- 

The radial part of this wave function V'ooHPo) can be expressed by a series of 5-functions: 

CVoHE^^fc^, (12) 

where pj(j = 1, . . . , N) are the grid points associated with the Gaussian nodes over the 
interval [—1, 1] with dimension equal N and Pn+i = Po- The coefficients C(j) are deter- 
mined from the solution of the linear algebraic equations system approximately equivalent 
to the Lippmann-Schwinger equation for two-neutron scattering [14J. 

In such a way we get the following expression for the amplitude of the nd charge- 
exchange process [6]: 

J = JpwiA + <J 1 S , 



Jpwia = - < LM L 1M S \IM D > x 



x {< lm 2 -m 3 \lM s >< m 1 m 2 ,Pi,Po + q/2\t° - t x |p,Po - q/2, mm 2 > x 

X ML (|p -q/2|)F^(p ^q/2)- 

- < ^rn' 2 -m 2 \lMs >< mim 3 ,pi,p - q/2\t° - t^p, p + q/2, mm' 2 > x 

x u L (\p + q/2\)Y^(p + q/2)} , (13) 



Ji So = V ; S m2 _ m3 < LM L 1M S \1M D >< -m"- - m' 2 \lM s > x (14) 

Z ■ 47T Z Z 

x J dp 'p ' 2 < m x m'z, Pi, p(, + q/2|t° - t^p, p' Q - q/2, mm" > VCfao) x 
x ML (|p[,-q/2|)r L ^(p^q/2) . 

Since the g,po "C p,Pi, we can neglect the g- and po- dependences of the high-energy 
np t-matrix in Eqs. (fT51) - (TT^l) . Besides, the integrand in Eq. tJTi]) is suppressed at high p' Q . 
It offers us the opportunity to consider np — > pn vertex as the free np scattering at angle 
6* = 7r in the center-of-mass system (c.m.s.) (or np —>■ pn at #i a b = 0). As well known, 
the NN t-matrix is described in collinear geometry by the three independent amplitudes: 

W = * ) = M® + [F(^) - B{E)\{°i&){°i?) + 5(S)(ffx<r a ), (15) 

where p* is the unit vector in the beam direction in the c.m.s. and E is the on-shell energy 
of two nucleons. Note, that this description is correct only for on-energy shell t-matrix in 
the center-of-mass system. However, we need the off-energy shell t-matrix in the frame, 
which in our kinematics corresponds to the laboratory one. In order to relate two these 
t-matrices, we use some results of the relativistic potential theory pjj [16] . This procedure 
has been presented in details in [5]. Here we give only the final formula: 

< mim2,ppo|t|ppo, mm" > = NN'F < mi\D^(u, p)|//i > x (16) 

x < m' 2 \D^(u,p' )\fi' 2 >< /ii/i2,p*|t cm |/i/i",p* > x 
x < n\D(u, p)|m >< /i"|L>(u, Po)|m" >, 

where D is a Wigner rotation operator in the spin space and u is a four-velocity. In 
our kinematical situation, where p = pi ^> p' , the Wigner rotation operators can be 
considered as the unity. In the other words, the spin structure of the tATAr-matrix in the 
reference frame is the same as the spin structure of the tjvTV-matrix in the c.m.s. 

The normalization factors N and N' are determined by the Jacobians for the transfor- 
mations between the reference frame and c.m.s. In our situation these factors are equal 
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to each other: 

-I 1/2 



N = N' 



m N + E p 



2m,N(mN + E p 



(17) 



4m N E p 

Other coefficient F is the kinematical factor connected with the transition from the on- 
energy shell i-matrix to off-energy shell one 



m N + E p 
2E p 



(18) 



Since we have neglected the pg dependence of the high-energy np t-matrix, Eq. (|14p can be 
integrated over p' taking into account Eqs.flE]), (1121) . In order to simplify final expressions, 
we consider here only S'-component of the deuteron wave function. This assumption 
does not to influence on the results, since we do not consider in the present paper any 
polarization observables, for which .D-wave contribution is very important. 

m N + E p 1 

x 



4:E p Vivr 
1/1 



{< ^m' 2 -m 3 \lM d >< m 1 m 2 ,p*\t° cm -t 1 cm \p*,mm l 2 > m (|Po - q/2|) - 
< ^2^ m 2|lM d >< mim 3 ,p*|t° m - tl m \p*,mm' 2 > u {\p + q/2|) + (19) 



+ {-l) l - m2 - m '*5m 2 , -ma < ~ m',\lM d > X 



1 „1 

— m - 

2 2 

N+1 Ci „ fa*+p* + q 2 /4" 



Here Qq(z) is the Legendre polynomials of the second kind. 

The cross section of the nd — > pnn reaction is defined in the standard manner: 

a = ( 2 7t) 4 ^ • 1 / d Pl dp 2 S(M d + E P -E 1 -E 2 ~ E 3 ) \J\ 2 , (20) 
p 6 J 



where p 3 = p — pi — p 2 and E 3 = \Jm? N + (p — Pl — p 2 ) 2 - The squared amplitude can 
be obtained straightforwardly from Eq.( jT9l) : 

\J? = ^ ( mN A E p Ep ) {( 252 (^) + i7l2 (^))[ 3M o(|Po-q/2|) + 3^(| P o + q/2|)- 



- 2u (|po-q/2|)u (|po + q/2|) + 
: m \ I'j'I 



+ |^£^( °? + ^ + ^ 4 W(fl| + (21) 



+ £ — Qo {^±A±U±\ [C ^(j) + C 001 (j)*] x 

; 'ii', V m J 
x M|po-q/2|) + u (|po + q/2|))] + 

+ 3A 2 (E)[u (\p + q/2|) - n (|p - q/2|)] 2 }. 

This expression contains both the spin-dependent part of the np amplitude, B and F 
terms, and the spin-independent one, A amplitude. However, A amplitude is multiplied 
by the difference of the two deuteron wave functions, which depend on the practically 
undistinguished arguments. The Eq. pT]) can be significantly simplified by the assumption 
of q = 0: 

\J\ 2 - ± (^^) 2 (2fi 2 (E) + F 2 (E))| Mo ( P o)+E 1 u Q { Pj )C^(j)\ .(22) 

In this expression the term corresponding to the spin-independent part of the np — > pn 
amplitude has vanished. Due to that we get the factorization of the squared nd — > 
pnn amplitude in the two parts. One of them depends on the deuteron and two slow- 
neutron wave functions. Other term corresponds to the spin- dependent component of the 
elementary np — > pn cross section. This result is very important, since it offers us the 
opportunity to extract the spin-dependent part of the neutron-proton charge-exchange 
cross section from more complicated reaction with the deuteron participation. But it 
should be noted, that such factorization is possible due to special kinematical conditions, 
when the transfer momentum is small in respect to the beam one. Moreover, we have used 
some model to take account of the two slow neutron final-state interaction. Therefore, 
the obtained result is mo del- dependent, what does not allow us to correctly extract the 
spin dependent part of the np cross section. 

3. Results 

In order to relate our calculation with the existing experimental data [lj we consider 
the ratio of the nd charge-exchange differential cross section to the free np scattering 
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differential cross section at 9^ = 0: 

^ da(nd — > pnn) ^da{np — > pn) 
dpidQ dQ 

The recent energy-dependent phase shift analysis data (PSA) [I7J have been used for 
the determination of the np amplitudes, which are needed to define both the nd and np 
charge-exchange differential cross sections. 

This ratio for the initial neutron kinetic energy T n = 1 GeV is presented in Fig. 2 as 
a function of the final-proton momentum pi. The dashed and solid curves correspond 
to the PWIA and P WIA+FSI calculations, respectively. One can see, the behaviours of 
these curves are significantly distinguished. The solid line has a very sharp peak, when 
the momentum p\ is close to the beam momentum p, or transfer momentum q is close 
to zero, while we do not observe any peak for the dashed line. This peak indicates the 
FSI contribution to the nd differential cross section. In this region the value of the R 
ratio varies in 10 times, when transfer momentum changes on a few MeV/c. Since any 
experiment has the limited momentum resolution, we consider also the R ratio integrated 
over pi in some region: 

Jp-Ap Jp-Ap dpidil dil 

Here we introduce a new variable Ap, which is a small difference between the initial 
neutron momentum and outgoing proton one. The integration limits change from p — Ap 
up to maximal value of p\ equal to p. 

The integrated R ratio is shown in Fig 3. in dependence on Ap at the neutron kinetic 
energy T n = 1 GeV. One can see, the PWIA curve is close to the PWIA+FSI one, 
when Ap increases. As it follows from Fig. 2, the FSI contribution at first increases and 
then decreases the nd differential cross section in respect to the PWIA predictions. This 
influence of the final-state interaction on the nd differential cross section has an effect on 
the behaviour of the integrated R ratio. In fact, the difference between the PWIA and full 
calculation results is about 30% for Ap equal to 10 MeV/c, about 15% for Ap equal to 20 
MeV/c and these lines are practically undistinguished, when Ap is equal to 60 MeV/c. 
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The energy dependence of the integrated R ratio is presented in Fig.4. The integration 
has been performed for Ap equal to 30 MeV/c. We investigate the energy region between 
800 and 1300 MeV. The upper limit is defined by the existing phase shift analysis data for 
np scattering. The dash-dotted line is obtained using the result of [H] with NN amplitudes 
taken from energy- dependent PSA [17]. The difference between result, obtained taking 
into account FSI, and PWIA result is about 10% for kinetic energy 800 MeV and few per 
cent for kinetic energy 1300 MeV. Thus, the contribution of the FSI decreases, when the 
kinetic energy is increasing. 

4. Conclusion 

In this paper the nd — > p(nn) reaction has been studied at the neutron kinetic energy 
T n = 0.8 — 1.3 GeV in kinematics, where transfer momentum is close to zero. It was shown 
that the expression for the nd — > pnn differential cross section is factorizable in the two 
parts. One of them is fully defined by the spin- dependent part of the elementary np — > pn 
cross section, while the other part depends on the deuteron and two slow neutrons wave 
functions. This factorization allows us to extract the spin-dependent part of the np charge- 
exchange squared amplitude, using the nd — > p{nn) reaction. But obtained result will 
be dependent on the model, which was applied for FSI description, and the choice of the 
deuteron wave function. This fact does not offer the opportunity to get the precise value 
of the spin-dependent part of the free np scattering amplitude by such method. However, 
it is possible to extract some useful information about np charge-exchange process (for 
example, sign, approximate value, etc.). 

The other important question, which has been studied, is the role of the final-state 
interaction of the two slow neutrons. We have considered the R ratio of the nd differential 
cross section to the elementary np — > pn one. It was shown, that the final-state interaction 
gives a considerable contribution into the three-fold nd — > pnn differential cross section, 
while the FSI influence on the integrated variables is small. Moreover, the contribution 
of the FSI into the integrated R ratio decreases with the increasing energy. 
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Figure captions 



Fig.l. Graphic representation of the amplitude of the nd — > pnn reaction. 

Fig. 2. R ratio vs. the fast proton momentum pi at T n — 1 GeV. 

Fig. 3. Integrated R ratio as a function of Ap at T n = 1 GeV. 

Fig. 4. Energy dependence of the integrated R ratio. 
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Figure 1: Graphic representation of the amplitude of the nd — > prm reaction. 
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P 1 , GeV/c 

Figure 2: R ratio vs. the fast proton momentum p 1 at T n — 1 GeV. 
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Figure 3: Integrated R ratio as a function of Ap at T n — 1 GeV. 



17 




18 



